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We present Floquet-Keldysh non-equilibrium dynamical mean field calculations (DMFT) with an 
iterative perturbative solver (IPT) for high pump-power induced Stark states yielding a closing gap 
within semiconductor bulk and nano-cavities. Our model predicts unusual broadening of bands in 
systems like ZnO random lasers. This can be explained as an exciton-polariton coupling within the 
siirgle crystalline pillar meetiirg the Mie resonance condition. Extraordinary enhancement of the 
lifetime of electronic states within the gap can lead to stable lasing modes and gain narrowing. 


The research of random and polariton lasing has fas¬ 
cinated the community ever since their discovery [ll-[l^. 
Both phenomena of coherent emission originating from 
semiconductor micro-resonators have been interpreted as 
to be fundamentally different. We show in this letter 
that on the basis of ultrafast non-equilibrium excitation 
the groundbreaking implication is possible to interpret 
random lasing as coupled cavity polariton lasers. On 
the other hand it would also allow to deduce that semi¬ 
conductor cavity polariton lasers develop laser thresholds 
which is the signature of a phase transition. 

In this letter we consider the semiconductor bulk as 
insulator, which is true in the zero temperature limit 
(T = 0), where all lattice influences are frozen out. The 
gap of 3.3 ek of ZnO corresponds to a Coulomb repul¬ 
sion U of about U/D = 4.0, which is deeply in the 
Mott insulating regime. D is the half bandwidth. First, 
the non-equilibrium bandstructure (i) of a classically 
strong field pumped bulk ZnO semiconductor is inves¬ 
tigated by Keldysh-Floquet non-equilibrium dynamical 
mean field theory (DMFT) and solved using second or¬ 
der non-equilibrium iterative perturbation theory (IPT). 
Second we explore topological effects (ii) induced in the 
non-equilibrium bandstructure when it is coupled to the 
cavity mode of a geometrical resonator (Fig. [2]). By com¬ 
parison of the life-times of both non-equilibrium solutions 
we identify long-lived resonances, exciton-polaritons, in 
the cavity coupled solutions which match in the frequency 
range experimentally determined spectra of both polari¬ 
ton lasers in ZnO nano-pillars as well as ZnO random 
lasing experiments. 

The Hubbard-type Hamiltonian Eq. m describes the 
electronic dynamics responsible for the generation of mi¬ 
croscopic gain and thus leading in further consequence to 
electronic inversion. It is generalized for electron-photon 
interaction which is also solved by applying the Floquet- 
Keldysh formalism using dynamical mean field theory 
(DMFT) and an iterative perturbative solver (IPT) (see 



FIG. 1: Insulator to metal transition, (a) Split energy band. 
The local Coulomb interaction U determines the width of the 
gap symmetrical to the Fermi energy Ef- (b) Electrons in the 
crystal lattice structure (black). Periodic driving yields an ad¬ 
ditional hopping contribution T{t) and local renormalization 
of the lattice energy efc = eo + E{t). The colors represent 
three time steps of driving. 


Fig. P. The classically field driven cavity-coupled 
Hamiltonian in position space reads 
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The first term ^ denotes the onsite pro¬ 
cesses which is completed by ^ <t 4 -^G.-cr de¬ 

voted to the Coulomb interaction U between electrons 
with opposite spins. The third term 
denotes the hopping processes with the amplitude t be¬ 
tween nearest neighbored sites. The last three terms 
are devoted to photons, including the external field- 
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FIG. 2: (a) The semiconductor behaves in this regime as an 
insulator: In the non-equilibrium DMFT-scheme optical exci¬ 
tation yields electronic hopping processes. They are mapped 
to the interaction with the single site on the background of 
the surrounding lattice (bath), (b) Integration over all sites 
leads to an effective theory including non-equilibrium excita¬ 
tion. This approach here is self-consistent because the bath 
consists of these very same single sites. The driven electronic 
system is additionally coupled to a cavity-resonance, the Mie 
resonance, (c) Experimentally relevant setup, a disordered 
waveguide of ZnO nanostructures embedded in a Si waver. 


electron coupling, the kinetic energy of the cavity pho¬ 
tons and the cavity photon-electron interaction. In 

id- Eo COsi^LT) J2<ij> {4,aCj,a “ 

renormalization of the hopping processes due to classical 
interaction of the dipole moment d with an external time 
dependent electric field of the amplitude Eq , the pumping 
laser, r denotes its time dependency and the corre¬ 
sponding laser frequency. The last two terms denote the 
cavity mode hujoCi'a at the resonance frequency coq as 
well as the coupling of this resonance to the electronic 
density of states gJ2k a 4 coupling 

strength g. The creator (annihilator) of a cavity pho¬ 
ton is denoted a (a’^') and of an electron it is c(c’f). The 
subscript k denotes the site, i,j denote next neighbored 
sites (Fig. [1]) within the originally assumed single-band. 
The first four terms of Hamiltonian Eq. (HD describe (i) 
externally pumped bulk semiconductor. 

The explicit time dependence of the external field yields 
Green’s functions that depend on two separate time ar¬ 
guments. A double Fourier transform from time- to fre¬ 
quency space leads to relative and center-of-mass fre¬ 
quency |13l. Il 4 | 
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{m,n) label the Floquet modes and (a,/3) address the 
branch of the Keldysh contour (±) where the respec¬ 
tive time argument resides. Floquet modes describe the 
Fourier-transformation of a periodic potential in time as 
the principal structure of bands in frequency space. How¬ 
ever the physical consequence which is noteworthy is a 
quantized absorption and emission of energy quanta hVlL 
of and to the original classical field. For non-interacting 
electrons the analytical solution for Gmn{k,uj) is solving 
the equation of motion. Photo-induced hopping leads to 
the retarded Green’s function for this sub-system 


G^^ik,u:) = J2 
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where ik represents the dispersion relation induced by the 
external field and is to be distinguished from e the lattice 
dispersion. The J„ are the cylindrical Bessel functions of 
integer order, Aq = d ■ Eq and is the laser frequency. 
The physical Green’s function is given according to 


GSQik,uj) = Y,G^„ik,uj). (4) 

m,n 


The Hamiltonian is solved through single-site Floquet- 
Keldysh DMFT (IPT) [3, [l^. This requires strict pe¬ 
riodicity of the lattice and homogeneous driving. That 
condition implies Bloch states and is characterized by 
the wave vector k within the laser-band-electron Green’s 
function G^Q(k,uj). Due to these requirements and the 
additional cavity in our system the DMFT selfconsis¬ 
tency relation, the centerpiece of DMFT, assumes the 
form of a matrix equation of dimension 2 x 2 in Keldysh 
space of the nonequilibrium Green’s functions and of 
dimension n x n in Floquet space. As DMFT solver 
we generalized the iterated perturbation theory (IPT) 
to Keldysh-Floquet form as well. The resulting numer¬ 
ics, although challenging, proves to be efficient and sta¬ 
ble also for all values of the Goulomb interaction U. 
Whereas we consider in 1^ a renormalized lattice po¬ 
tential, we take into account here a coupling of the mi¬ 
croscopic dipole moment to the external field amplitude 
^ Il 4 | . This implies the quantum-mechanical character 


of the dipole operator d = id J2<ij> {4,aCj,a “ 
in the fourth term of Eq. O- It is fundamentally differ¬ 
ent from the generic kinetic hopping of the third term. 
The coupling d ■ AoCOs(DiT) under the assumption of 


the Coulomb gauge E{t) = —-^A{t), which reads in 
Fourier space E{V,l) = i^L ■ A{Ql), generates the fac¬ 
tor Hi that cancels the I/Hi term in the renormalized 
cylindrical Bessel function in Eq. (7) of ref. [l^. Con¬ 
sequentially we consider the first n = 10 Floquet modes 
only which has been proven to be sufficient here. 

The physical consequences are fascinating since in 
the non-equilibrium externally pumped state the con- 
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FIG. 3: LDOS for ZnO bulk semiconductor. In equilib¬ 
rium the electronic gap is S.SSeV, the microscopic dipole mo¬ 
ment equals |d| = 2.30826288 x the lattice con¬ 

stant a = 0, 325 nm. The pump power is increased from 0 to 
5.0MW/cm^. A significant shift of states into the semicon¬ 
ductor gap is observed. This energy region has been experi¬ 
mentally identified to be the gain area for polaritonic states. 
Corresponding life times can be found in Fig. [I] 


FIG. 4: Inverse lifetime 1/r corresponding to the LDOS (Fig. 
0 for increasing pump power, (a) The peak centered at 
1.69eF indicates small lifetimes, fast decay of mid-gap states 
(red emission), (b) Finite lifetime of non-equilibrium states 
near the band edge indicate stable states which can couple to 
geometric (Mie) resonances and form cavity-polaritons. With 
rising pump power the lifetime drops but remains larger than 
the lifetime of nearby excitonic states in the band. 


ventional bulk semiconductor ZnO exhibits a counter¬ 
intuitive behavior. The parameters are chosen such that 
the resulting band split equals the ZnO semiconductor 
band gap of 3.3 eV in the non-pumped case. Through 
classical external excitation a multi-band insulator de¬ 
velops, where the quantized exchange manifests itself 
in a multitude of sub-bands. The original bandgap of 
3.3 eF is subsequently occupied with dressed states or 
excitons [l^. The complex permittivity is defined as 
e = e [l — = e' — ie" where e' = CrCo is the real part 

of the permittivity. The optical conductivity determines 
the imaginary part we" = a/cue 

Qp2x2 p 
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and characterizes the microscopic gain spectrum (see 
refs. 11, Isl). 

The optical cavity and the Fano-coupling ITj to the elec¬ 
tronic subsystem of the bulk semiconductor is described 
(ii) in the last term of the original Hamiltonian, Eq. O- 
It yields an additional contribution to the self-energy 
E in the denominator of the retarded Green’s function 
Eq.(0. The form of E significantly triggers the numer¬ 
ical efficiency of the selfconsistency consideration. The 
strength g of the coupling between optically excited elec¬ 
trons ^ and the cavity resonance hujoa^a has to 

be considered g < 1 in order to guarantee convergence 
of the perturbation series up to second or higher order. 


Our results for the non-equilibrium bandstructure (i) of 
the classically excited ZnO bulk (Fig. [3]) show for in¬ 
creasing pump strengths (colored graphs) spectral weight 
inside the semiconductor band gap whereas for zero exci¬ 
tation (black curve) the original semi-circular bandstruc¬ 
ture for non-interacting systems (band-width assumed to 
be 1.7 eV) splits through the Coulomb interaction U in 
upper and lower band. The tails at the outer upper and 
lower bandedges are typical for results derived by the 
DMFT (IPT) solver. It is found that the additional Flo- 
quet bands emerge with increasing pump intensity. Fur¬ 
ther the development of moderate sub gaps towards the 
near-gap band edges at —0.67 eV and 4.05 eV (compared 
to the bandedge of the valence band) is seen whereas the 
band edges themselves show a sharp step towards the 
novel spectral weight inside the gap. Interestingly the 
relation of the sub-gaps and the simultaneously found 
features in the results for the non-equilibrium lifetimes 
1/ImE of the excited states (see Fig. 0]) is non-trivial. 
We find for the same excitation wavelength (A = 710nm) 
and a raising intensity the increase of the non-equilibrium 
lifetime at the Fermi edge (see Fig. 0^). Also 1/t is larger 
in the gap near the band edges which corresponds with 
a higher local density of electronic states (LDOS), i.e. a 
rising mobility. A plateauing of lifetimes is observed near 
the edge, whereas deep in the band at 4.8eE we derive 
for increasing intensity a crossover. Here 1/r develops 
a peak towards |ifg| = ‘i.iMW/cm? which flips into a 
sharp dipped feature (red dash-dots) and arrives at a 
lifetime almost as high as in the non-excited case (black 
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FIG. 5: Polaritonic lifetimes. Pump frequency A = 710 nm, 
cavities resonance A = 505.34 nm, refractive index n = 
1.97182, pump intensity 2AMW/cm~^. The difference to 
pure bulk material is evident in novel, sharp features (marked) 
at 0..0.5eP and 2.5..3.3el7. These energies correspond to 
spectrally separated lasing peaks in random laser systems [T^ . 
Tuning the coupling strength of bulk matter and cavity yields 
small effects. Inset: Optical DC-conductivity of the semicon¬ 
ductor Mie resonator for 7 = 0.2 representing processes. 

line). This behavior is rather unexpected however it is 
confirmed by a consistency check of the Floquet sum. 

As results for the cavities (Fig. [S]) we derive profound 
changes compared to pure bulk which are only in parts 
expected. The nano-pillar shall support a single whisper¬ 
ing gallery mode, (a) creates (annihilates) a photon as 
the energy eigen-state at ?iuJo= 2,45el7 (= 505.34nm). 
This is the in-plane Mie resonance, which is related to 
the continuum of dressed electronic states (Fig. 2]) and 
so the non-equilibrium complex refractive index of ZnO. 
The momentum of the cavity photon is small compared to 
the electron’s momentum and (/photon — 0 is set whenever 
we consider the electronic subsystem. We compare the 
lifetime for varying coupling strength (g/t) = 0.2, 0.4, 0.9 
at zero temperature and at half filling, yielding a suppres¬ 
sion of the spectral function around the Fermi-level (half 
width of the gap) where electrons are transferred to the 
high (low) energy tails of the spectral function. Our re¬ 
sults are derived for a pumping intensity of 2AMW/cm^ 
which is typical for solid state random lasers. In agree¬ 
ment with results for the bulk we also find the charac¬ 
teristic feature around the Fermi edge. Also the dip at 
4.8 eV persists in the cavity system. The physical inter¬ 
pretation of the non-equilibrium lifetimes including the 
Mie resonance in (Fig. S]) leads to ultrafast electro-optical 
Kerr effects, non-equilibrium |xp processes with are of 
the order of 0.75ps up to 1.25ps. This is a long-living 
cavity polariton which couples to the non-equilibrium 
AC Stark effect of the pure excited bulk. The peaks 
in the lifetime of the semiconductor cavity around 2.8eV 
and 3.2eV are occurring comparably far from the cavi¬ 


ties resonance frequency and so they are a proof of the 
importance of the Fano resonance for electronic spec¬ 
trum in general and they express the splitting in several 
relevant polariton branches which define the lasing fre¬ 
quencies. This is a property of the independent single 
scatterer, which as for the random lasing experiments, 
only profits from randomness by the increased impinging 
photonic intensity. The optical DC-conductivity (inset 
FiglU) of the semiconductor Mie resonator for 7 = 0.2 
in the non-equilibrium, i.e. |xp processes in the non¬ 
equilibrium, shows a strong frequency-dependency and 
remarkable features which have not been reported be¬ 
fore to the best of our knowledge. The dips may lead to 
significant gain narrowing. 

Our results suggest, that solid state random lasers can 
be interpreted as coupled cavity polariton lasers where 
the Mie cavity plays a crucial role in the generation 
of topological features in the microscopic electronic 
bandstructure of the semiconductor. We theoretically 
demonstrated polaritonic gain in ZnO as non-equilibrium 
effect in the Hubbard Model. The polaritonic states 
are comprised of long-lived excitonic AC Stark states in 
the original ZnO gap and the Mie resonances of ZnO 
nano-pillars. This implies that the ZnO laser plasma 
is actually a polariton condensate which may display a 
threshold behavior. To distinguish between the concept 
of either cavity polaritons or excitons and the particle 
hole plasma and additionally the influence of the mode, 
i.e. the quality and the coupling of the cavity, a temper¬ 
ature dependent DMFT in the the non-equilibrium will 
be useful. This novel theoretical approach leads to a 
more fundamental understanding of ultrafast topological 
effects in semiconductor nano-cavities. 
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